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Abstract. We investigate photon emission from a moving particle in an expanding 
universe. This process is analogous to the radiation from an accelerated charge in 
the classical electromagnetic theory. Using the framework of quantum field theory 
in curved spacetime, we demonstrate that the Wentzel-Kramers-Brillouin (WKB) 
approximation leads to the Larmor formula for the rate of the radiation energy from a 
moving charge in an expanding universe. Using exactly solvable models in a radiation- 
dominated universe and in a Milne universe, we examine the validity of the WKB 
formula. It is shown that the quantum effect suppresses the radiation energy in 
comparison with the WKB formula. 



1. Introduction 

One of the notable feature of quantum fields in curved spacetime is that quantum 
processes prohibited in the Minkowski spacetime is allowed [H [2]. For example, the 
emission of a photon from a moving massive charged particle occurs in an expanding 
universe, though such a process is prohibited by the energy momentum conservation 
in the Minkowski spacetime due to the Lorentz invariance. This subject was studied 
by several authors: Pioneering work was done by Buchbinder and Tsaregorodtsev [3] 
and by Lotze[4j. These authors investigated the transition probability of the process by 
applying QED to a radiation-dominated universe. Then Futamase et al. and Hotta et 
al. considered the case of a simple background spacetime with a sudden transition of the 
scale factor [3, [6]. These previous works, however, focused on the transition probability 
of the process. While, in the present paper, we calculate the radiation energy emitted 
through the process. Our point of view is as follows: The motion of a massive charge in 
an expanding or contracting universe can be regarded as an accelerated motion, because 
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the physical momentum of the particle decreases (increased) as the universe expands 
(contracts). Then, the photon emission process can be regarded as the well-known 
classical radiation process from an accelerated charge [7]. The present paper aims to 
clarify the correspondence between the classical and quantum approaches to photon 
emission from a moving charge in expanding universe. 

This paper is organized as follows: In section 2, we review the scalar QED model in 
the Friedmann-Robertson- Walker universe. Then, we show that the classical radiation 
formula, which corresponds to the Larmor formula for the rate of the radiation energy 
from an accelerated charge, is derived under the WKB approximation in the framework 
of quantum field theory. In section 3, we analyze the scalar QED model without WKB 
approximation in two different universe models. We first consider a universe which 
undergoes a bounce and is radiation-dominated at asymptotic past and future infinity. 
The advantage of this model is that the equation of motion of a free complex scalar field 
is exactly solvable. Then, the radiation energy through the photon emission process 
is numerically computed. The result is compared with the corresponding result based 
on the WKB formula, and the condition for the validity of the WKB formula is found. 
Second, we consider the case of a bounced Milne universe, to verify the robustness of the 
result. The last section is devoted to the summary and conclusions. Throughout this 
paper, we use the unit light velocity equals 1. We follow the convention (— , +, +, +). 

2. Derivation of the radiation formula with the WKB approximation 

In what follows, we focus on the spatially fiat Friedmann-Robertson- Walker spacetime, 
whose line element is expressed as 

rfs^ = a{viY[-drf + dx^] = aijiYri^^dx^dx" , (1) 

where r] is the conformal time, and 0(77) is the scale factor. We consider the scalar QED 
Lagrangian conformally coupled to the curvature. 



S 



J d'^Xy/^ 



2 +econfi? F^'^Fm- 



V ^ / 4/io 

where F^,^ = V^Ai^ — V^A^ is the field strength, and /io is the magnetic permeability 
of vacuum. In this paper, we explicitly include the Planck constant h. Introducing the 
conformally rescaled field 0, 



we may rewrite the Lagrangian as 



S= d^x 



h J \ h J h 4/io 

where /^,^ = d^A^ — dyA^. Thus the system is mathematically equivalent to the scalar 
QED in the Minkowski spacetime with the time- variable mass ma{T]). 
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Figure 1. Fcynman Diagram of tlic photon emission process 



We follow a general prescription for interacting fields, based on the interaction 
picture approach (see e.g., [H 12]). We focus on the radiation energy emitted through 
the process described by the Feynman diagram in Fig. 1. The S-matrix corresponding 
to the diagram is 

^ = T / ^ (05^0^ - d'x. (2) 

Here the field operator in the right-hand side, A'^ and (p obey the free field equations. 
The quantized free fields are expressed as follows. For the photon field, we have 



M^) = E E ^/?^^ ^ \ ^ + Hermite Conjugate \ , (3) 



r=l,2 fc 



where uj{k) = k, e^^^ is the polarization vector, and ajj'' and a^j^^'' are the annihilation 
and creation operators. They satisfies the commutation relations. 



"fc ) "fc' 



The vacuum state of the photon field is 

a^k^\0)^^'^ = 0, for any k, r. 
The complex scalar field is quantized as 



"fc > "fc' 



0. 



v^eI 



where hq and ¥ satisfy 



6, 



bq, hq, 



=0, 



(4) 



(5) 



(6) 



and Cq and satisfy the same commutation relations. The mode function ^qirf) satisfies 
the equation of motion 



+ 



with the normalization condition 



df] 



dr] 



I. 



(7) 



(8) 



The vacuum state of the complex scalar field is 

b^\0f^ = c^\0f^ = 0, for any q. (9) 

It may be noted that because of the time- variation of the mass ma{T]) the definition 
of the vacuum state can be ambiguous. In this section, however, we assume that there 
exists a natural stable vacuum state, which is the case when the WKB approximation 
is valid. 

With the WKB approximation, the mode function (Pq{ri) is given as 



where 



2^],(r/) 



exp 



—I 



h 



We can write the condition that the WKB formula is valid, as (see e.g., [T]) 
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(10) 
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(12) 



where the dot means the differentiation with respect to rj. 

Using the WKB mode functions, we evaluate the transition amphtude of the process 
described by Fig. 1, 

Transition Amplitude = (fl^li) , (13) 

where the initial state and the final state are 

|i) = 6tjo)«|0)(^), (14) 

\i)=KaVW'W'\ (15) 

respectively. 

The radiation energy is 
(27r)3 



E 



L3 



^ ^ /ico'(A;)|Transition Amplitude|" 



r=l,2 k q 



which gives 



Bo J (27r)3 



(16) 



(17) 



where qf = qi — k, and Eq is the permittivity of vacuum, which is related to /io as 
EofiQ = 1/c^ = 1. We set qi ■ k = qikcosO, and rewrite (ITTl) as 



2eo 



d^k 



X 



drj 



(27r)3"* 
1 



qUi-cos'e 



exp 



ikr] + i / dr]'^lq^. {rj') - i / drj'Qq^ {rj') 



(18) 



where 



(19) 



where the first and last equahties follow from the definition of the momentum of the 
charged particle, 



p = hq. 

We follow the prescription in ref. [S], 



dq 



Iq2 + ^Hty 



and 



where we used qf = qi — k and 

k{=\k\) <^qf{=\qf\), and qi{= \q,\) . 



Then, we have 



E 



2eo ^0 (27r)2 J^i 



d cos oil — cos 9 



X 



dv^T^ exp 



ik / dr]' [1 



Pi cos 9 ' 



We can choose pi to be in proportion to 2;-axis. Therefore, we may write 
, dz n dr] 

Pi=p = -r^ p =-r^ 

dr dr 
where r is the parameter chosen as 
dx^ dx" 



-(p) +Pz = -m-^a^T]) 



dr dr 

Because we may write Pi/p^ = dz/drj, we have 



Jdv 



Pi 



p^iv) 



exp 



dr]'k 1 



Pi cos 9 ' 



p^{'ri') 

= J dzexp [ik(r](z) — cos 9z 

Instead of z or rj, introducing the variable ^ defined by 
= T] — cos 9z{ri), 

we have 



E = — , 

2eo Jo 



r°° dkk"^ /"i / r, \ r dz 



(27r 



dz 



(20) 



(21) 
p=Pi 



(22) 
(23) 



(24) 



(25) 



(26) 



(27) 
(28) 

(29) 
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Furthermore, the partial integration gives 



. (30) 

where we have assumed that the boundary term has no contribution. The neglect of 
the boundary term is justified if \dz/d^\ is finite at all times. Since we have 

d^ l-COS^i' ^ /p2 + ^2^2 



where z = dz/dr], we see that \dz/dC^ \ is bounded at all times. 
The integration with respect to k yields 

where the integration is extended to — oo < A; < oo, and divided by the factor 2 [8]. 
From the definition of ^, we have 

d^ dj^ z z 

— = {l- zcosO), J7^ = T. : H\3- (22) 

drj dt, (1 — -2;cosy)'^ 

Then, Eq. (15^ is rephrased as 

2 \ ■■2 

E=-^[ dcosen - cos^ e) [ dri- — . (34) 

Finally, the integration with respect to cos 6 yields 

This result is the same as the Larmor formula in the case when the particle moves along 
a straight line [7]. In the present case, Eq. (135!) is rewritten as 

Note that this is the energy in the conformally rescaled spacetime, which is not the 
physical energy. From this expression, however, we can read the physical radiation 
energy as S = E/a, and the physical rate of the radiation energy per unit time as 

d£ ^ ]_dE ^ eVpf,y,H^ 
dt a? drj QireQim? ' 
where t is the cosmic time, Pphys{= Pi/d) is the physical momentum and H is the Hubble 
expansion rate. This result is consistent with that in ref. [5], which is obtained from 
consideration of the classical electromagnetic radiation formula of a moving charge in 
an expanding universe. 

Finally in this section, let us summarize the necessary condition for reproducing 
the radiation formula in the classical electromagnetic theory. We started with the WKB 
formula of the mode function, for which Eq. ( fT2l) is needed. In addition, we assumed 
Eq. (l23l) . Although this condition is independent of the necessary condition for the 
WKB approximation, Eq. (|T2l) . we assumed it because this additional assumption is 
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Figure 2. The scale factor of the bounced radiation dominant universe as a function 

•n- 

necessary to recover the conventional picture for the classical radiation from a charged 
massive particle, in which the massive field should behave like a particle and the photon 
field should behave like a wave. 

3. Quantum radiation formula 

In this section, we calculate the radiated energy without the WKB approximation. We 
start with reviewing the formalism. Now we return to the S-matrix expression ([2]). A 
characteristic feature of the quantum field theory in curved spacetime is that vacuum 
state may not be stable. Namely, the vacuum state may change, and the particle creation 
phenomenon can occur. This vacuum effect is not taken into account in the analysis 
based on the WKB approximation in the previous section. In this section, we adopt the 
formalism for evaluating the S-matrix element, taking the vacuum effect into account. 

For definiteness, we assume that the interaction is switched off in the asymptotic 
past infinity [i] = —oo) and future infinity [rj = +oo), where the different vacuum states 
1 0)1^-* and |0)q^^, respectively, are defined for the free field. Then, similar to Eq. ([5]), we 
may write the quantized field as 



or as 




where 6™, and 6^°^, c™^ are the annihilation and the creation operators, respectively, 
with respect to the in-vacuum aX rj = — oo, while 6°*^*, c°"* and 6°"*^, c"*^*^ are those 
with respect to the out- vacuum at ?7 = +oo. The in-vacuum and out- vacuum states are 
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expressed as 



and 



= = 0, (40) 



^riO)iut = ClO)iut = 0, (41) 

for any q. These annihilation and the creation operators satisfy the same commutation 
relations as Eq. ([6]). 

In Eqs. fl38|) and fl39|) . v^g (//) and v^q"*(?7) are the mode functions with respect to 
the vacuum states \0)[? and |0)^^^, respectively. They are related by the Bogoliubov 
transformation, 

<(r^) = a,<Hr^) + (42) 
where and (3q satisfy the normalization condition 

i^^r - = 1- (43) 

The creation and annihilation operators are related as 

K = «;^r - (44) 

ci_V = a,c-*t_;3^5-\ (45) 
Using the above relations, we have 

vTi^K'^ + vTivr^^ = ^^^C^ + ^^^c% (46) 



and 



out/ \ m^'^*(n) 

VaivK + ^r(^)cLV = + :^^c°-t. (47) 



As for the photon field, because of its conformal invariance, the vacuum state is invariant 
in time, |0)(^) = |0);^) = |0)« . 

The transition amplitude in the lowest order of e is evaluated in the similar form 
as Eq. f[T^ . but with 

|i)in = ^'j:^|0)^'^|0)!?, (48) 

|f)out = ^r^«?^|0)^'^|0)i?, (49) 

In the computation of the transition amplitude, it is necessary to regularize the 
divergence arising from the vacuum-to-vacuum amplitude. In the flat background, this 
can be done unambiguously by taking the normal-order product of operators. However, 
in a curved spacetime, there arises ambiguity because the in-state annihilation/creation 
operators are different from the out-state annihilation/creation operators. In fact, there 
will be particle creation from vacuum and the vacuum-to-vacuum amplitude will not 
longer be unity any more. 

To deal with this situation properly, we consider the generalized normal product of 
operators, which is defined as the form where the operators are expressed only in terms 
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of the in-state annihilation operator and the out-state creation operators, and all the 
out-state creation operator are placed to the left of all the in-state annihilation operators 
[9]. This is adopted in ref. [3]. The nice properties of this generalized normal ordering are 
that it is symmetrically defined with respect to the in-state and out-state operators, and 
that the vacuum-to-vacuum amplitude is normalized to unity. In particular, this latter 
property means that it minimizes the effect of particle creation from vacuum. Since what 
we are interested in is not the vacuum particle creation but the transition amplitude 
for a massive particle to radiate under the electromagnetic interaction, we adopt the 
generalized normal ordering to regularize the divergence. We note that what would 
be actually observed will be inevitably contaminated by the effect of particle creation. 
However, whether there is a way to separate out this vacuum effect observationally is 
beyond the scope of the present paper. 
Following ref. [3], we define 



out(f|iV[5]|i), 



Transition amplitude = — j— — —. (50) 

Here |0)j^ = |0)^^^|0)J^^ and \0)^^^ = |0)^^^ lO)^^^, respectively, N[S] denotes the 
generalized normal product of the S-matrix, and out(0|0);j^ is the vacuum to vacuum 
amplitude. 

We then find the formula for the radiation energy as 



{2TT)'^eoJo J-l Ittq.PI". 



9/1 



(51) 



3.1. Example I 

In this subsection, we consider a time-symmetric bounce universe which asymptotically 
approaches a contracting and expanding radiation-dominated universe (see Fig. 2). The 
scale factor is given in terms of the conformal time rj by [lOl [11] 



a{ri) = Y f)2j^2 _|_ g2 (— c)o < 7] < oo), (52) 

which recovers a radiation-dominated Friedmann universe in the asymptotic regions, 
a{t) (X t^/"^ {r] ±oo). In this background spacetime, the WKB formula fl35|) gives 

Equation of motion ([7]) is reduced to the Weber's differential equation 
where 

"=2^ ±rn.,lk - = ±(l±')vW7S',. (55) 



1 



F 



0.1 



1 1 










■ 


■ 


exact 




■ ■ ■ ■ approximate 
1 1 







0.1 1 

p. /me 

Figure 3. F as function of Pi/me with fixed as hg/me^ ~ O.Of , 0.1, 1, from top to 
bottom. The solid curve is exact, Eq. (|66|) . while the dotted curve is approximate, 
Eq. (pj). 




Figure 4. F as function of fig/rae'^ with fixed as Pi/me = 0.01 (dashed curve) and 
1 (solid curve). 
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Therefore, the mode functions are constructed as 

1 



e-"^/^AA-i/2(-(l - i),Jmg/hr]), 



with 



X 



2mg/h 

Using the mathematical formula for the parabolic cylinder function 



DJz) 



we easily find the Bogoliubov coefficients 



Or 



,7r(-A+j/2)/2 



(3q = e 



n{-X-i/2) 



" r(-zA+l/2) 

Using the mathematical formula for the parabolic cylinder function [13 



(56) 
(57) 

(58) 
(59) 
(60) 



/oo 
-oo 

^exp 
u(j,a, 1 + i{a — P), 



4 4 ^ ^^4^ ^ a/2 



X 



(61) 



we have 



\/2mg 



+ — A - A + + z A - A In 

4 4 ^ ^ Amg/h ^ ' 



X 



f/(a + ^,l + (A-A), 



2'' ' 2mQin 

where t/(a, c, z) is the second order confluent hypergeometric function 
the notation 



A 
A 



2mQ/?i 2 gh \ 

~ 2^ 
o Pi 



2mg/h 



1 me^ 



Pt 



2 gh 



m?e^ me me^ 



Pi \ 










{p^-nkf^ 




k cos 9 





2mg/h 
(62) 
and we used 

(63) 



Then, the radiation energy E is represented by 

Pi T^Q 



(64) 
(65) 
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where 



Pi hg \ 12 f°° 



- em e^m 



TC Jo 



dkk^ 



dcosO- 



1 - cos^ 9 



,7r(A-2A) 



-2ttX 



X 



)(1 
1 



c/(a + j.i + .(A-A).-.M: 



(66) 



and we have defined k = {e/g)k. Note that the function F describes the deviation from 
the WKB formula. 

Now let us consider the classical limit by taking the limit h ^ 0. We use the 
mathematical formula 

22(i-)/2K,_i(2v/i) 



lim U(a, c, z/a) , , , 

where Ky{z) is the modified Bessel function. Then, we have 



(67) 



Pi hg \ 24 /-"o 



dkk^ 



dx{l - x2)e"(^-^) 



X 



i(X-X) 



k\l + 



Pt 



with 



1 hg ^2 



A - A = —kx ^ „ 
me I me^ 



k'. 



(6^ 



(69) 



We have evaluated the function F numerically. The cases of Pi/me = 0.01 and 1 are 
shown in Fig. |H In all cases we analyzed, we found F is a decreasing function of hg/me^. 
In the limit — > 0, 

Pi 



A-A 



-kx = s . 



Then, we have 



me 



24 



lim F = — 



dkk' 



dx{l - x^)e''' 



KiAk\\^ 



Pi 



rp?e^ 



(70) 



(71) 



Thus F is the function of only Pi/me in this limit. By numerical analysis of the right- 
hand side of Eq. (I7T!) , we find 

(72) 



lim F 



irrespectively of pi/ me. We therefore infer that the integral of Eq. ( iTTj) gives 1. Although 
we have not yet succeeded to show it in general, we can show that, for the case 
Pi/me <C 1, Eq. (17T|) yields 

24 /•! 



lim F ~ — / dx(l - x"^) / dkk^Ko(k? = 1, 

h^O TT^ 7-1 Jo 

by using the integral formula for the modified Bessel function |14j . 

x"^ K ^{ax)Ky{hx)dx 



(73) 



r(l-A) 
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^•1 — A + /i — i^N^^^l — A — /i — 1/ 

X 




2 ^ ^ 2 

xF(i^^i^,i^^^;i-A;i-y. (74) 

This demonstrates that the exact formula in the hmit of ^ agrees with the WKB 
approximate formula. Then, the decrease of F from 1 comes from the term in proportion 
to ^, hence the suppression is understood as the quantum effect. 

Let us summarize the result. The quantum radiation formula agrees with the WKB 
formula under the condition 

(75) 

^ (76) 
me 2 me^ 

We can show that the former condition flTSl) is derived from the condition for the WKB 
approximation, Eq. (fT2l) . while the latter condition (!76l) is equivalent to Pi ^ hk, i.e., 
Eq. (l23l) . The former condition is satisfied when the Compton wavelength of the charged 
particle is shorter than the Hubble horizon length around the bounce regime defined by 
a ~ e ~ ^ (see below), where the classical radiation rate becomes maximum. Figure [3] 
plots F as a function of Pi/me with fixed as hg/me^ = 1, 0.1 and 0.01. Figure H] plots 
F as a function of hg/me^ with fixed as Pi/me = 0.01 and 1. These figures show F = 1 
for Pi/me ^ 1 and hg/me^ <^ 1, and the suppression F < 1 for the other region. 

Finally in this section, we mention the physical meaning of these parameters. 
Around the bounce regime a ~ e ~ grj, we can write 

Pi Pphys 



J77,£ fYi bounce 7716 bounce 

where Pphys = Pi/^ is the physical momentum, Xq = h/m is the Compton wavelength, 
and -ff|bounce = is the Hubble expansion rate, respectively. Thus pi/me is the 
relativistic factor, and hg/me^ can be regarded as the the ratio of the Hubble horizon 
length to the Compton wavelength of the charged particle, around the bouncing regime. 

3.2. Example II 

Here we consider a Friedmann spacetime with the scale factor 
/ +i+e (t>0) ^ r ee+^ (r^ > 0) 

[ -i+e (t<0) [ ee-^ (r/ < 0) ' ^ ^ 

where e(> 0) is a small constant, i = t/to and f/ = ri/rjQ. This mimics a Milne-like 
bounce universe but with a flat spatial geometry. 

The solution of the Klein-Gordon equation is written with the Bessel function, and 
the positive frequency mode function is [15] 

H%imi+t + e)) (t>0) 
H^l\m{-t + e)) (t<0) 



iPg = NgX'{ (if , (79) 
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where 

Ng = e-^-/'^ (80) 

from the normahzation condition, and we have defined m = riQin/h and q = rjoq. 
The in- and out-vacuum mode functions are therefore given by 

a,H%imii+e)) + f3,Hil\m{t + e)) {t > 0) 

Hff{m{-t + e))) {t<0) ^ ^ 



X J H%{m{t + e))) it>0) 
' > agHil\m{-t + e)) + PgH%{m{-t + e)) {t < 0) 



(82) 



where the coefficients are 



«.-^e---(<(me))-, 

= -^e-'^^-(i/«(m6)//i|(m6))', 



nme 



(3, = ^e-'^«-(//i^^(me)i/i%(me))' = (83) 

and the prime means Hi^z)' = dH{z)/dz. 
We may write 

f (me)-*^(m(t + e))*^ (f > 0), 

\ {mey''{m{-t + e))"^^ (t < 0), 



then, we have 



(84) 



roo 



/oo 
dr]~ 
-oo 

+ r]oN,^N,^{mef' H dzz-^'-'H%^{z) {a,,H^^^z) + P,,H%^{z)) . 



(85) 



In the hmit z ^ 1, we have 

/^W(^) = ^e^^+^'^/2-^./^ 

Thus in the hmit, me = mrjoe/h ^ 1, we may approximate 

= ^ -^e2^'^^ (86) 
P,. = ^ 0. (87) 
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and we have 

2 

rf7^e^^Vf,Vr = - ze^--^- _ \ (88) 
After the integration with respect to cos 6', the total energy is 

E = r S-^^ = -J^. (89) 



This completely agrees with what is obtained from the integration of fl36|) with 
Note that this result is obtained under the condition me = mrjQe/h ^ 1, which may be 
expressed as 

me = mrjoe/h = Xc'H-\^^^^^ » 1, (90) 

where Ac = h/m is the Compton wavelength and -ff|bounce = l/(e'7o) is the Hubble 
expansion rate at the bounce regime a ~ e, i.e, {\ri\ ^ rjo). Thus the WKB formula is 
valid as long as the Compton wavelength is shorter than the Hubble horizon length. 



4. Summary and Conclusions 

In the present paper, we investigated photon emission from a moving massive charge in 
an expanding universe. We considered the scalar QED model for simplicity, and focused 
on the energy radiated by the process. First we showed how the Larmor formula for the 
rate of the radiation of energy in the classical electromagnetic theory can be reproduced 
under the WKB approximation in the framework of the quantum field theory in curved 
spacetime. 

We also investigated the limits of the validity of the WKB formula, by deriving 
the radiation formula in a bouncing universe in which the mode functions are exactly 
solvable. The result using the exact mode function shows the suppression of the radiation 
energy compared with the WKB formula. The suppression depends on the ratio of the 
Compton wavelength Ac of the charged particle to Hubble length H^^. Namely, the 
larger the ratio \c/H~^ is, the stronger the suppression becomes. In the limit the 
Compton wavelength is small compared with the Hubble length, the radiation formula 
is found to agree with the WKB formula. Since this limit is equivalent to the limit ^ — 0, 
the suppression we found is a genuine quantum effect in an expanding (or contracting) 
universe, which is due to the finiteness of the Hubble length. Whether the quantum 
effect on the radiation from a accelerated charge always leads suppression or not is an 
interesting question. This would be understood by analyzing higher order terms of the 
WKB approximation. We will return to this point in future. 
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